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PHYSICAL PROPERTIES OF THE 0.94-DAY PERIOD TRANSITING PLANETARY SYSTEM WASP-18* 
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ABSTRACT 

We present high-precision photometry of five consecutive transits of WASP-18, an extrasolar planetary sys- 
tem with one of the shortest orbital periods known. Through the use of telescope defocussing we achieve a 
photometric precision of 0.47-0.83 mmag per observation over complete transit events. The data are analysed 
using the JKTEBOP code and three different sets of stellar evolutionary models. We find the mass and radius of 
the planet to be M b = 10.43 ± 0.30 ± 0.24 M Jup and R b = 1 . 165 ± 0.055 ± 0.014R Jup (statistical and systematic 
errors) respectively. The systematic errors in the orbital separation and the stellar and planetary masses, aris- 
ing from the use of theoretical predictions, are of a similar size to the statistical errors and set a limit on our 
understanding of the WASP-18 system. We point out that seven of the nine known massive transiting planets 
(Mb > 3 Mj up ) have eccentric orbits, whereas significant orbital eccentricity has been detected for only four of 
the 46 less massive planets. This may indicate that there are two different populations of transiting planets, but 
could also be explained by observational biases. Further radial velocity observations of low-mass planets will 
make it possible to choose between these two scenarios. 
Subject headings: 



1. INTRODUCTION 

The recent disco very of the transit ing extrasolar planetary 
system WASP-18 dHellieret all 120091 hereafter H09) lights 
the way towards understanding the tidal interactions between 
giant planets and their parent stars. WASP-18 b is one of 
the shortest-period (P or b = 0.94 d) and most massive (Mb = 
10Mj up ) extrasolar planets known. These properties make 
it an unparallelled indicator o f the tidal dissipation param- 
eters (Goldreich & Soter 1966) for the star and the planet 
(I Jackson et al.ll2009l) . A value si milar to that observed for So- 
lar system bodies (Q ~ 10 4 -10 6 : |Peale1l 19991) would cause the 
orbital period of WASP-18 to decrease at a sufficient rate for 
the effect to be observable within ten years (H09). 

In this work we present high-precision follow-up photom- 
etry of WASP-18, obtained using telescope-defocussing tech- 
niques dSouthworth et al.l 12009 a b) which give a scatter of 
only 0.47 to 0.83 mmag per observation. These are analysed 
to yield improved physical properties of the WASP-18 sys- 
tem, with careful attention paid to statistical and systematic 
errors. The quality of the light curve is a critical factor in 
measurement of the physical properties of transiting planets 

* BASED ON DATA COLLECTED BY MINDSTEP WITH THE DANISH 
1 .54 M TELESCOPE AT THE ESO LA SILLA OBSERVATORY 
' Royal Society University Research Fellow 



dSoufhworfhll2009l) . In the case of WASP-18 the systematic 
errors arising from the use of theoretical stellar models are 
also important, and are a limiting factor in the understanding 
of this system. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed five consecutive transits of WASP-18 on the 
nights of 2009 September 7—11, using the 1.54m Danish 
Telescope at ESO La Silla with the focal-reducing imager 
DFOSC. The plate scale of this setup is 0.39" pixel" 1 . The 
CCD was windowed down to an area of 12' x 8' in order to 
decrease the readout time to 51s. This window was chosen 
to include WASP-18 (HD 10069, spectral type F6 V, V = 9.30, 
B— V = 0.44) and a good comparison star (HD 10179, spectral 
type F3 V, V = 9.65, B-V = 0.45). 

A Johnson V filter and exposure times of 80 s were used (in- 
stead of our usual Cousins R and 120 s), to obtain lower count 
rates from the target and comparison stars. We defocussed the 
telescope to a point spread function (PSF) diameter of ^90 
pixels (35"), to limit the peak counts to 45 000 per pixel for 
WASP-18 and also lower the flat-fielding noise. This focus 
setting was used for all observations and the pointing of the 
telescope was maintained using autoguiding. An observing 
log is given in TableQ] 

We took several images with the telescope properly fo- 
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Table 1 

Log of the observations presented in this work. N b s is the number of observations. 'Moon' and 'Distance' are the fractional illumination of the Moon at its 

distance from WASP- 1 8 at the midpoint of the transit. 



Date Start time (UT) End time (UT) N \, a Exposure time (s) Filter Airmass Moon Distance (°) Scatter (mmag) 



2009 09 07 06:41 10:05 93 80.0 V 1.05^1.27 0.860 61.9 0.83 

2009 09 08 05:00 10:04 140 80.0 V 1.15 -> 1.04 1.28 0.784 67.1 0.68 

2009 09 09 03:55 10:00 169 80.0 V 1.31 -> 1.04 -> 1.27 0.696 73.3 0.51 

2009 09 10 03:25 08:29 141 80.0 V 1.41 -> 1.04 -> 1.10 0.597 80.2 0.56 

2009 09 11 04:07 07:12 86 80.0 V 1.25^ 1.04 0.492 87.6 0.47 



Table 2 

Photometric observations of WASP-18. This table is included to show the 
form of the data; the complete dataset can be found in the electronic version 
of this work. 



Midpoint of 


Relative R 


Error in R 


observation (BJD) 


magnitude 


magnitude 


55082.783033 


-0.000200 


0.000788 


55082.784491 


-0.001340 


0.000790 


55082.786030 


-0.000495 


0.000786 


55082.787570 


0.000825 


0.000786 


55082.789109 


0.000989 


0.000804 


55082.790649 


0.001469 


0.000880 


55082.792188 


-0.000431 


0.000965 


55082.793727 


0.000483 


0.001159 


55082.795267 


0.001142 


0.001067 


55082.796806 


-0.001431 


0.001156 



cussed in order to check that there were no nearby stars con- 
taminating the PSF of WASP-18. There are five detectable 
objects nearby, but all are at least 100 pixels distant. The 
brightest is 6.75 mag fainter than WASP-18 and the other four 
are more than 9.5 mag fainter. We conclude that the PSF of 
WASP-18 is not contaminated by any star detectable with our 
equipment. 

Data reduction was pe rformed in the same way as in 
ISouthworth et al.1 (l2009albh . In short, we run a reduction 
pipeline written in I DL 1 which use s an implementation of the 
DAOPHOT package (Stet sor3[l987l) to perform aperture pho- 
tometry. For each night's data the apertures were placed in- 
teractively and their positions fixed for each CCD image. We 
found that the most precise photometry was obtained using an 
object aperture of radius 60 pixels and a sky annulus of radius 
85-110 pixels. The choice of aperture sizes (within reason) 
excludes flux from the nearby faint stars, and has a negligible 
effect on the resulting photometry. 

We calculated differential magnitudes for WASP-18 using 
an ensemble of four comparison stars, of which HD 10179 
is by far the brightest. The light curve for each night was 
normalised to zero differential magnitude by fitting a straight 
line to the observations taken outside transit, whilst simulta- 
neously optimising the weights of the four comparison stars. 
We have applied bias and flat-field corrections to the images, 
but find that this does not have a significant effect on the pho- 
tometry. The individual light curves are shown in Fig. [2 and 
the full 629 datapoints are given in Table|2] The scatter in the 
final light curves varies from 0.47 to 0.83 mmag per point, and 
is higher for data taken when the moon was bright (TableQ]). 

3. LIGHT CURVE ANALYSIS 

1 The acronym IDL stands for Interactive Data Language and is a 
trademark of ITT Visual Information Solutions. For further details see 

http : //www . ittvis . com/Product Services / IDL . aspx 
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Figure 1. Light curves of WASP-18 from the five nights of observations. 
For each night the error bars have been scaled to give a reduced x OI 
x\ = 1-0. 



The light c urves of WASP-18 were analysed with the JK- 
TEBOP 2 code (ISouthworth et al.ll20Q4allbh. u sing the approach 
discussed in detail by ISouthworthl (120081) . jktebop ap- 
proximates the two components of WASP-18 using biaxial 

2 JKTEBOP is written in FORTRAN77 and the source code is available at 

[http : / /www, astro . keele . ac . uk/~ jkt /| 
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Figure 2. Phased light curve of WASP-18 and the best fit found using JKTEBOP and the quadratic LD law. The residuals of the fit are offset from zero to appear 
at the base of the figure. 



spheroids whose shapes are governed by the mass ratio; we 
adopt the value 0.008 but find that large changes in this num- 
ber have a negligible effect on the results. 

WASP-18 b has a slightly eccentric orbit, which should be 
taken into account as it has a small effect on the shape of 
the transit. However, for transit light curves the ef fect is in 
gener al too subtle to include it as a fitted parameter (Kipping 
2008). We could fix eccentricity, <?, and periastron longitude, 
oj, to the values obtained from the velocity variation of the 
parent star (H09), but this would neglect their measurement 
uncertainties. We have therefore modified JKTEBOP to allow 
the inclusion of e and u as fitted parameters constrained by the 
known values and uncertainties. In practise we use ecosw = 
0.0008 ± 0.0014 and e sin oj = -0.0093 ± 0.0030, as these two 
parameters are only weakly correlated with each other. In the 
current case the uncertainties in the ecosw and esinw values 
have a very minor effect on our results. 

We have incorporated the H09 photometry of WASP-18 in 
order to obtain the most precise ephemeris. This was done by 
including the reference transit epoch from H09 as an observed 
quantity, an d then fitting for transit epoch and orbital period 
directly (see South worth et al.ll2007l) . We chose a new transit 
epoch which is close to the midpoint of our own observations, 
so is essentially based on just the data presented in this work. 
The final eclipse ephemeris is 

T = BJD 2455084.792931(88) + 0.94145181(44) x E 

where Tq is the transit midpoint, E is the number of orbital 
cycles after the reference epoch, and quantities in parenthe- 
ses denote the uncertainty in the final digit of the preceding 
number. 

The limb darkening (LD) of WASP 1 8 A was accounted fo r 
using five different parametric laws (see ISouthworthl 2008). 
Theoretical LD coefficients were obtained by bilinear inter- 
polation, to the known effective tempera ture (T B ff) and surfac e 
gravity of the parent star, in the tables oflVan Hammd d 19931) . 
lai^dlM 120043) and lClaret & Hauschildj (120031 We ob- 
tained solutions with the LD coefficients fixed at the theoreti- 
cal values, with the linear coefficient fitted for and the nonlin- 



ear coefficient fixed (and optionally perturbed by ±0.05 on a 
flat distribution in the error analyses), and with both LD coef- 
ficients included as fitted parameters. The full set of solutions 
is given in Table[3] 

The uncertainties of the light curve par ameters were as- 
sessed using Monte Carlo simulations (Sout hworth et alj 
l2004cl l2005bl) . The importance of red noise was checked 
using a residual-permutation approach, and found to be mi- 
nor. The solutions with LD coefficients fixed to theoretical 
values are poorer than those where one or both LD coeffi- 
cients are fitted parameters. We adopt the mean of the so- 
lutions with non-linear LD and both LD coefficients fitted, 
as these are the most internally consistent. The final uncer- 
tainties come from Monte Carlo solutions but include contri- 
butions from the residual-permutation analyses and the (mi- 
nor) variation between the solutions with different LD laws. 
We find the fractional radius 3 of the star and planet to be 
r A = 0.2795 ± 0.0084 and r b = 0.0272 ± 0.0012, respectively, 
and the orbital inclination to be i = 85.0° ± 2.1°. 

4. THE PHYSICAL PROPERTIES OF WASP-18 

The physical properties of a transiting planetary system 
cannot in general be calculated purely from observed quan- 
tities. The most common way to overcome this difficulty is to 
impose predictions from theoretical stellar evolutionary mod- 
els onto the parent star. W e have used tabulated prediction s 
fro m three sources: Claret dClaretll2004bl 120051 120061 [ 2007), 
Y 2 dDemaraue et al]l2004l) and Cambridge liPols et alJl 19981 
Eldridge & Tout 2004). This allows the assessment of the sys- 
tematic errors caused by using stellar theory. 

We began with the parameters measured from the light 
curve and the observed velocity amplitude of the parent star, 
K A = 18 18.3 ±8. 0ms -1 (H09). These were augmented by an 
estimate of the velocity amplitude of the planet, to calcu- 
late preliminary physical properties of the system. We then in- 
terpolated within one of the grids of theoretical predictions to 

3 Fractional radius is the radius of a component of a binary system ex- 
pressed as a fraction of the orbital semimajor axis. The utility of this quantity 
is that it is measureable from light curve data alone. 
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Table 3 

Parameters of the JKTEBOP best fits of our V-band light curve of WASP- 18, using different approaches to LD. For each part of the table the upper quantities are 
fitted parameters and the lower quantities are derived parameters. ta and rj> are the fractional radii of the star and planet, respectively, and k = r^/rj^. i is the 
orbital inclination, and ha and va are the linear and non-linear LD coefficients, respectively. P is the orbital period and Tq, the reference epoch of minimum 

light, is given as BJD - 2455000.0. 





Linear LD law 


Quadratic LD law 


Square-root LD law 


Logarithmic LD law 


Cubic LD law 


All LD coefficients fixed 




0.3094 ± 0.0075 


0.3034 ± 0.0074 


0.3060 ± 0.0077 


0.3069 ± 0.0075 


0.3176 ±0.0068 


k 


0.09704 ± 0.00075 


0.09636 ±0.00057 


0.09682 ± 0.00069 


0.09712 ±0.00062 


0.09912 ±0.00048 


i (deg.) 


84.68 ± 1.33 


85.96 ±1.70 


85.31 ± 1.51 


85.05 ± 1.35 


83.05 ±0.90 


"A 


0.60 fixed 


0.40 fixed 


0.20 fixed 


0.70 fixed 


0.40 fixed 


VA 




0.30 fixed 


0.60 fixed 


0.25 fixed 


0.15 fixed 


P 


0.94145182 ± 0.00000040 


0.94145181 ±0.00000041 


0.94145182 ± 0.00000044 


0.94145181 ±0.00000045 


0.94145181 ±0.00000043 


Tq 


84.792935 ± 0.000090 


84.792929 ± 0.000085 


84.792932 ± 0.000084 


84.792931 ±0.000084 


84.792924 ± 0.000085 


>'A 


0.2821 ±0.0067 


0.2767 ± 0.0066 


0.2790 ± 0.0069 


0.2798 ± 0.0067 


0.2890 ± 0.0060 


rb 


0.02737 ± 0.00085 


0.02667 ± 0.00078 


0.02701 ± 0.00085 


0.02717 ±0.00081 


0.02864 ± 0.00072 


a (mmag) 


0.6406 


0.6358 


0.6356 


0.6355 


0.6408 




1.0194 


1.0048 


1.0031 


1.0034 


1.0247 


Fitting for the linear LD coefficient and fixing the nonlinear LD coefficient 




0.3106 ±0.0070 


0.3053 ± 0.0082 


0.3073 ± 0.0080 


0.3066 ± 0.0076 


0.3071 ±0.0079 


k 


0.09813 ± 0.00070 


0.09672 ± 0.00087 


0.09720 ± 0.00082 


0.09702 ±0.00081 


0.09726 ± 0.00086 


i (deg.) 


84.2 ± 1.1 


85.5 ±1.8 


85.0± 1.5 


85.1 ±1.5 


85.0± 1.5 


"A 


0.527 ± 0.023 


0.384 ±0.029 


0.1 80 ±0.024 


0.705 ± 0.025 


0.495 ±0.025 


VA 




0.30 fixed 


0.60 fixed 


0.25 fixed 


0.15 fixed 


P 


0.94145182 ±0.00000043 


0.94145181 ±0.00000041 


0.94145182 ± 0.00000042 


0.94145181 ±0.00000042 


0.94145182 ±0.00000039 


To 


84.792931 ±0.000080 


84.792930 ± 0.000089 


84.792931 ±0.000086 


84.792931 ±0.000081 


84.792932 ± 0.000084 


I'A 


0.2829 ± 0.0062 


0.2784 ± 0.0072 


0.2801 ±0.0071 


0.2795 ± 0.0068 


0.2799 ± 0.0070 


rb 


0.02776 ± 0.00078 


0.02692 ± 0.00094 


0.02722 ± 0.00090 


0.02711 ±0.00086 


0.02722 ± 0.00088 


a (mmag) 


0.6357 


0.6356 


0.6352 


0.6355 


0.6350 


x\ 


1.0048 


1.0061 


1.0037 


1.0049 


1.0031 


Fitting for the linear LD coefficient and perturbing the nonlinear LD coefficient 




0.3106 ±0.0066 


0.3053 ±0.0084 


0.3073 ± 0.0079 


0.3066 ±0.0081 


0.3071 ±0.0079 


k 


0.09813 ±0.00069 


0.09672 ± 0.00092 


0.09720 ±0.00081 


0.09702 ± 0.00086 


0.09726 ± 0.00083 


i (deg.) 


84.2 ± 1.0 


85.5 ±1.8 


85.0± 1.5 


85.1 ±1.6 


85.0±1.5 


"A 


0.527 ± 0.022 


0.384 ±0.039 


0.1 80 ±0.042 


0.705 ± 0.053 


0.495 ±0.028 


VA 




0.30 perturbed 


0.60 perturbed 


0.25 perturbed 


0.15 perturbed 


P 


0.94145182 ± 0.00000042 


0.94145181 ±0.00000043 


0.94145182 ± 0.00000043 


0.94145181 ±0.00000043 


0.94145182 ± 0.00000043 


Tq 


84.792931 ±0.000086 


84.792930 ± 0.000085 


84.792931 ±0.000080 


84.792931 ±0.000084 


84.792932 ± 0.000084 


ta 


0.2829 ± 0.0059 


0.2784 ± 0.0075 


0.2801 ±0.0070 


0.2795 ± 0.0072 


0.2799 ± 0.0070 


rb 


0.02776 ± 0.00076 


0.02692 ± 0.00095 


0.02722 ± 0.00086 


0.02711 ±0.00090 


0.02722 ± 0.00091 


a (mmag) 


0.6357 


0.6356 


0.6352 


0.6355 


0.6350 


X V 


1.0048 


1.0061 


1.0037 


1.0049 


1.0031 


Fitting for both LD coefficients 


TA + r-b 


0.3106 ±0.0070 


0.3064 ± 0.0086 


0.3068 ± 0.0083 


0.3063 ± 0.0083 


0.3071 ±0.0081 


k 


0.09813 ± 0.00070 


0.09740 ±0.001 17 


0.09721 ±0.00144 


0.09737 ±0.00128 


0.09721 ±0.00140 


i (deg.) 


84.2 ± 1.1 


85.0 ±1.7 


85.0± 1.9 


85.0 ±1.7 


85.0± 1.8 


"A 


0.527 ±0.022 


0.473 ± 0.087 


0.213 ±0.433 


0.617±0.175 


0.492 ± 0.045 


VA 




0.11 ±0.18 


0.54 ±0.77 


0.12 ±0.24 


0.16 ±0.20 


P 


0.94145182 ±0.00000041 


0.94145181 ±0.00000044 


0.94145181 ±0.00000043 


0.94145182 ± 0.00000044 


0.94145181 ±0.00000041 


To 


84.792931 ±0.000086 


84.792932 ± 0.000089 


84.792931 ±0.000084 


84.792932 ± 0.000088 


84.792931 ±0.000084 


>'A 


0.2829 ± 0.0062 


0.2792 ± 0.0076 


0.2796 ± 0.0072 


0.2792 ± 0.0073 


0.2799 ±0.0071 


rb 


0.02776 ± 0.00078 


0.02719 ±0.00103 


0.02718 ±0.00108 


0.02718 ±0.00103 


0.02721 ±0.00105 


a (mmag) 


0.6357 


0.6353 


0.6352 


0.6354 


0.6350 


X v 


1.0048 


1.0055 


1.0053 


1.0057 


1.0047 



find the expected radius and T e ff of the star for the preliminary 
mass and the measured metal abundance ( [||] = 0.00 ± 0.09; 
H09). Kb was then iterative ly refined to minimise the differ- 
ence between the model-predicted radius and T e s, and the cal- 
culated radius and measured T e « (6400 ± 100 K; H09). This 
was done for a range of ages for the star, and the best overall 
fit retained as the optimal solution. Finally, the above process 
was repeated whilst varying every input parameter by its un- 
certainty to build up a complete err or budget for each output 
parameter (So uthworth et al.ll2005al) . A detailed description 



of this process can be found in Southworth (200H). 

Table|4] shows the results of these analyses. The physical 
properties calculated using the Claret and Y 2 sets of stellar 
models are in excellent agreement, but those from the Cam- 
bridge models are slightly discrepant. This causes systematic 
errors of 4% in the stellar mass and 2% in the planetary mass, 
both of similar size to the corresponding statistical errors. The 
quality of our results is therefore limited by our theoretical 
understanding of the parent star. Our final results are in good 
agreement with those of H09 (Table|4]i, but incorporate a more 
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Table 4 

Physical properties for WASP- 18, derived using the predictions of different sets of stellar evolutionary models. For quantities with two errorbars, the 

uncertainties have been split into statistical and systematic errors, respectively. 



Cambridge models 



Y 2 models 



Claret models 



Final result (this work) 



Hellier et al. (2009) 



Orbital separation (AU) 

Stellar mass (Mq) 

Stellar radius (R©) 

Stellar log g [cgs] 

Stellar density ( pq ) 

Planetary mass (Mj U p) 

Planetary radius (Rjup) 
Planet surface gravity (ms" ! ) 

Planetary density (PJup) 

Stellar age (Gyr) 



0.02022 ± 0.00022 

1.235 ±0.039 
1.215 ±0.048 
4.361 ±0.022 
0.689 ±0.062 
10. 18 ±0.22 
1.151 ±0.052 
191 ±17 
6.68 ±0.89 
0.0-0.6 



0.02043 ± 0.00028 
1.274 ±0.052 
1.227 ±0.043 
4.365 ±0.026 
0.689 ±0.062 
10.39 ±0.30 
1.163 ±0.055 
191 ±17 
6.61 ±0.89 
0.0-2.1 



0.02047 ± 0.00027 

1.281 ±0.050 
1.230 ±0.042 
4.366 ±0.026 
0.689 ±0.062 
10.43 ±0.28 
1.165 ±0.054 
191 ±17 
6.60 ±0.88 
4+' 2 

-04 



0.02047 ± 0.00028 ± 0.00025 
1.281 ± 0.052 ±0.046 
1.230 ± 0.045 ±0.015 
4.366 ± 0.026 ±0.005 
0.689 ± 0.062 ±0.000 
10.43 ± 0.30 ±0.24 
1.165 ± 0.055 ±0.014 
191 ± 17 

6.60 ± 0.90 ±0.08 
0.0 - 2.0 



0.02026 ± 0.00068 
1.25 ±0.13 
1 216+ 0067 
4 367 +0.028 

707+ 056 

'-0.096 

10.30 ±0.69 

1 106 +0 .072 
i.iuo „ 054 

194+12 

7 7^+0.78 

0.5-1.5 



comprehensive set of uncertainties. 

These results give the equilibrium temperature of the planet 
to be one of the highest for the known planets: 

/1-A\ 1/4 
T eq = (2392 ±63) f —j K 

where A is the Bond albedo and F is the heat redistribution 
factor. This equilibrium temperature, and the closeness to its 
parent star, make WASP-18b a good target for the detection 
of thermal emission and reflected light. 

5. CONCLUSIONS 

We have presented high-quality observations of five con- 
secutive transits by the newly-discovered planet WASP- 18 b, 
which has one of the shortest orbital periods of all known tran- 
siting extrasolar planetary systems (TEPs). Our defocussed- 
photometry approach yielded scatters of between 0.47 and 
0.83 mmag per point in the final light curves. These data 
were analysed using the JKTEBOP code, which was modified 
to include the spectroscopically derived orbital eccentricity in 
a statistically correct way. The light curve parameters were 
then combined with the predictions of theoretical stellar evo- 
lutionary models to determine the physical properties of the 
planet and its host star. 

A significant source of uncertainty in our results stems from 
the use of theoretical models to constrain the physical proper- 
ties of the star. Further uncertainty comes from observed T e ff 
and [^] , for which improved values are warranted. However, 
the systematic error from the use of stellar theory is an im- 
portant uncertainty in the masses of the star and planet. This 
is due to our fundamentally incomplete understanding of the 
structure and evolution of low-mass stars. As with many othe r 
transiting systems (e.g. WASP-4; ISouthworth et al.l l2009b), 
our understanding of the planet is limited by our lack of un- 
derstanding of the parent star. 

We confirm and refine the physical properties of WASP 18 
found by H09. WASP- 18 b is a very massive planet in an ex- 
tremely short-period and eccentric orbit, which is a clear indi- 
cator that the tidal effects in planetary systems are weaker than 
expected (see H09). Long-term follow-up studies of WASP- 
18 will add progressively stricter constraints on the orbital de- 
cay of the planet and thus the strength of these tidal effects. 

We now split the full sample of known (i.e. published) 
TEPS into two classes according to planetary mass. The 
mass distribution of transiting planets shows a dearth of ob- 
jects with masses in the interval 2.0-3.1 Mj up . There are nine 



planets more massive than this and 46 less massive. Seven 
of the nine high -mass TEPs h ave eccentric orbits (HAT-P-2, 
Bakos et al.l2007t H P 17156, lBarbieri et al.l2007t HP 80606, 
Laughlin et ail2 009; WASP-10, Christi anTt"al.l2009t WASP- 
14. iJoshi et al.l l2009t WASP-18: XO-3. Uohns-Krull et al l 

2008) , and the existing radial velocity observations of the 
remaining two cannot rule out an e ccentr icity of 0.03 or 
lower (CoRoT-Ex o-2, lAlonsoetalJ 12001 OGLE-TR-L9, 
Snellen et al. 2009). By comp arison, only four of the 46 low- 
mass TEPs have a significant dLucv & Sweenevll 1971b orbital 
eccentricity measurement. 

These numbers imply that the more massive TEPs are a dif- 
ferent population to the less massive ones; Fisher's exact test 
(Fishen [l922l) returns a probability lower than 10~ 5 of the null 
hypothesis (although this does not account for our freedom to 
choose the dividing line between the two classes). This in- 
dicates that the two types of TEPs have a different internal 
structure, formation mechanism, or evolution, a suggestion 
which is supported by observations of misalign ment between 
the sp in and orbital axes of M > 3Mj up TEPs (Johnso n et alJ 

2009) . 

There is, however, a bias at work here. The more massive 
TEPs cause a larger radial velocity signal in their parent star 
(Mb oc K 3 ), so a given set of radial velo city measurements ca n 
detect smaller eccentricities (see also IShen & Turner! 120081) . 
The eccentricity of the WASP-18 system is in fact below the 
detection limit of existing observations of most TEPs. We 
therefore advocate the acquisition of additional velocity data 
for the known low-mass TEPs, in order to equalise the ec- 
centricity detection limits between the two classes of TEPs. 
These observations would allow acceptance or rejection of the 
hypothesis that more massive TEPs represent a fundamentally 
different planet population to their lower-mass brethren. 



The observations presented in this work will be made avail- 
able at the CPS ([http://cdsweb.u-strasbg.fr/ll 
and at |http : / /www . astro . keele . ac . uk/~> jkt/] 
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